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METHOD OF COMPUTING THE PITCH NAMES 
OF NOTES IN MIDI-LIKE MUSIC 
REPRESENTATIONS 

Field of the invention 

This invention relates to the problem of constructing a robust pitch spelling algorithm — 
that is, an algorithm that reliably computes the correct pitch names (e.g.. Cjf4. Bb5 etc.) 
of the notes in a passage of tonal music, when given only the onset-time, MIDI note 
number (The MIDI Manufacturers' Association, 1996, p. 10) and possibly the duration of 
each note in the passage. 

There are good practical and scientific reasons for attempting to develop a robust pitch 
spelling algorithm. First, until such an algorithm is devised, it will be impossible to con- 
struct a robust MIDI-to-notation transcription algorithm — that is, an algorithm that reli- 
ably computes a correctly notated score of a passage when given only a MIDI file of the pas- 
sage as input. Commercial music notation programs (e.g., Sibelius (www.sibelius.com). 
Coda Finale (www.codamusic.com) and Nightingale (www.ngale.com)) typically use 
MIDI-to-notation transcription algorithms to allow the user to generate a notated score 
from a MIDI file encoding a performance of the passage to be notated. However, the 
MIDI-to-notation transcription algorithms that are currently used in commercial music 
notation programs are crude and unreliable. Also, existing audio transcription systems 
generate not notated scores but MIDI-like representations as output (see, for example, 
Davy and Godsill, 2003; Plumbley et al, 2002; Walmsley, 2000). So if one wishes to 
produce a notated score from a digital audio recording, one typically needs a MIDI-to- 
not^tion transcription algorithm (incorporating a pitch spelling algorithm) in addition to 
an audio transcription system. 

Knowing the letter-names of the pitch events in 1 a passage is also indispensibie in music 
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information retrieval and musical pattern discovery (Meredith et at, 2002). For example, 
the occurrence of a motive on a different degree of a scale (e.g., C-D-E-C restated as 
E-F-G-E) might be perceptually significant even if the corresponding chromatic intervals 
in the patterns differ. Such matches can be found using fast, exact-matching algorithms if 
the pitch names of the notes are encoded, but exact-matching algorithms cannot be used 
to find such matches if the pitches are represented using just MIDI note numbers. If a 
robust pitch spelling algorithm existed, it could be used to compute the pitch names of 
the notes in the tens of thousands of MIDI files of works that are freely available online, 
allowing these files to be searched more effectively by a music information retrieval (MIR) 
system. 

In the vast majority of cases, the correct pitch name for a note in a passage of tonal 
music can be determined by considering the roles that the note is perceived to play in the 
perceived harmonic structure and voice-leading structure of the passage. For example, 
when played in isolation in an equal-tempered tuning system, the first soprano" note in 
Figure la would sound the same as the first soprano note in Figure lb. However, in 
Figure la, this note is spelt as a GJj4 because it is perceived to function as a leading note 
in A minor; whereas in Figure lb, the first soprano note is spelt as an A\>4 because it 
functions as a submediant in C minor. Similarly, the first alto note in Figure lb would 
sound the same as the first alto note in Figure lc in an equal-tempered tuning system. 
However, in Figure lb the first alto note is spelt as an Fl}4 because it functions in this 
context as a subdominant in C minor; whereas, in Figure lc, the first alto note functions 
as a leading note in Fjj minor so it is spelt as an E{j4. . 

Nevertheless, it is not always easy to determine the correct pitch name of a note 
by considering the harmonic structure and voice-leading structure of its context. For 
example, as Piston (1978, p. 390) observes, the tenor Eb4 in the thud and fourth bars of 
Figure 2 should be spelt as a D#4 if one perceives the harmonic progression here to be 
+II 2 - 1 as shown. But spelling the soprano Eb5 in the fourth bar as D#5 would result in 
a strange melodic line. 
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Such cases where it is difficult to determine the correct pitch name of a note in a tonal 
work are relatively rare — particularly in Western tonal music of the so-called 'common 
practice' period (roughly the 18th and 19th centuries). In the vast majority of cases, those 
who study and perform Western tonal music agree about how a note should be spelt in 
a given tonal context. This poses an interesting problem for cognitive science, namely: 
what are the cognitive processes involved when a musically trained individual determines 
the correct pitch name of a note in a passage of tonal music? 

A good way of trying to answer this question is to attempt to construct a robust pitch 

0 11 ^ ^ 1^~~:^U~>^ x-U^x. ;„ „ — i-U^ 4- :^ ^^^^^x-^^x. — ;±u ^/u^x. ;^ 1 — ^ „ — ^ U ^ , , x 

ojjciiiiig ai^uiniiiiii uxxcbu wjjct Cbocn ixx a vv<xy uuai ±0 owiiniouciiu vv i un vv xxct>u 1C5 Aiiwvvn auuuu 

the neurology, physiology, psychophysics and psychology of expert music perception and 
cognition. Such a pitch spelling algorithm could serve the valuable scientific purpose of 
furthering our understanding of the cognitive processes underlying an interesting intellec- 
tual ability exhibited by many individuals who study and perform Western tonal music. 

The vast majority of notes in authoritative published editions of scores of common 
practice tonal works are generally agreed to be spelt correctly by those who understand 
Western staff notation. Therefore a pitch spelling algorithm or computational model of 

pxu^xx t>pc;xxixi^ uclxx c v <xi ixohjCvx uujcouivci) L^y x ixxxxxixx^ 10 uii uwixa-i wuiao rx-xj ^jLiipadiu^ uxic 

pitch names it predicts with those of the corresponding notes in authoritative published 
editions of scores of the works. However, this can only be done accurately and quickly 
if one has^access to encodings of these authoritative scores in the form of computer files 
that can be compared automatically with the pitch spelling algorithm's output. 

> 

Related Art 

In this section, I compare the performance of three prior pitch spelling algorithms on a 
single test corpus of works containing 41544 notes and consisting of all 48 pieces in the 
first book of J. S. Bach's Das Wohltemperirte Klavter. The algorithms compared are those 
of Cambouropoulos (1996, 1998, 2000, 2001, 2002), Longuet-Higgins (1976, 1987, 1993) 
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and Temperley (1997, 2001). 

When these algorithms were run on the test corpus, Cambouropoulos's algorithm 
made 2599 mistakes, Longuet-Higgins's algorithm made 265 mistakes and Temperley's 
algorithm made 122 mistakes. 

The test corpus The test corpus used in the comparison consists of representations 
of the scores of all 24 Preludes and all 24 Fugues in the first book of J. S. Bach's Das 
Wohltemperirte Klavier encoded in what I call OPND format. 1 Each OPND representa- 
tion is a set of triples, (t, n, d), each triple giving the onset time, t, the pitch name, n and 
the duration, d of a single note (or sequence of tied notes) in the score 2 The onset time 
and duration of each note are expressed as integer multiples of the largest common divisor 
of all the notated onset times and note durations in the score. For example, Figure 3b 
lives the OPND representation of the score in Figure 3a. Note that within each pitch 
name in an element in an OPND -representation, each flat symbol is represented by an 
<F' character and each sharp symbol is represented by an <S' character (double-sharps are 
denoted by two 'S' characters, so, for example, Fx4 is denoted by "FSS4".) 

The test corpus was derived by automatic conversion from Hewlett's (1997) MuseData 
encodings of Bach's Das Wohltemperirte Klavier? The MuseData encodings contained 
minor errors which were corrected in the OPND representations. 

Also, Temperley's algorithm cannot deal with situations in which two or more notes 
with the same pitch begin at the same time. Thus, wherever two or more notes with the 
same pitch p began simultaneously at time t, all notes with pitch p and onset time t were 
removed from the OPND file except the one with the longest duration. This resulted in 
a test corpus containing 41544 notes. 4 

The "piano-roll" or MIDI-like input representations accepted by the algorithms com- 
1 OPND stands for "onset, pitch-name, duration" 

^ w^SS^^ * the mes used in the comparison ' but this voice information 18 

^ Available online at http://ww.musedata.org/encodings/bach/bg/keybd/ 
The complete test corpus is available online at http://www.titanmusic.com/data.html. 
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pared in this study were derived automatically from the OPND representations of the 
pieces in the test corpus. 

Longuet-Higgins's algorithm Pitch spelling is one of the tasks performed by Longuet- 
Higgins's (1976, 1987, 1993) music. p program. Longuet-Higgins has published the frill 
POP-11 source code of this program (Longuet-Higgins, 1987, pp. 120-126; Longuet- 
Higgins, 1993, pp. 486-492). Because of my lack of familiarity with POP-11, I decided to 
translate just the pitch spelling portion of music.p directly into Lisp in order to make it 
easier to perform the comparison. This was possible because the pitch spelling portion of 
music.p operates independently of the rhythmic part. 5 

The input to music.p must be in the form of a list of triples, (p, ton each triple 
giving the "keyboard position" p together with the onset time t on and the offset time £ Q ff 
in centiseconds of each note. The keyboard position p is simply an integer indicating the 
key that would have to be pressed on a normal piano keyboard in order to perform the 
note, with Ct|3 mapping onto 0, C|J3 and Db3 mapping onto 1, Cb4 mapping onto 12 and 
so on (i.e., p = MIDI NOTE NUMBER - 48). 

The algorithm then computes a value of "sharpness" q for each note in the input 
(Longuet-Higgins, 1987, p. 111). The sharpness of a pitch name indicates the position of 
the pitch name on the line of fifths (Temperley, 2001, p. 117) and is therefore essentially the 
same as Temperley's (2001, p. 118) concept of "tonal pitch class". In Longuet-Higgins's 
algorithm, if q x and q 2 are the sharpnesses of two notes then the interval between the 
notes is defined to have "degree" Sq = q<z — gi. If \6q\ < 6, the interval is defined to be 
"diatonic"; if \Sq\ > 6, it is defined to be "chromatic"; and if \Sq\ = 6, it is defined to be 
"diabolic" (Longuet-Higgins, 1987, p. 112). Longuet-Higgins's algorithm attempts to spell 
notes so that the degree between each note and the tonic at the point at which the note 
occurs is not chromatic (Longuet-Higgins, 1987, p. 113). The algorithm also incorporates 
various rules for dealing with chromatic passages (Longuet-Higgins, 1987, pp. 113-114). 

5 My Lisp implementation of the pitch spelling portion of music.p is available online at 
http: //www. titanimisic . com/soft ware .html. 
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Pitch spelling algorithms sometimes spell complete pieces in a different key from that 
in which they are notated in the original score. For example, Longuet-Higgins's algorithm 
spells the Prelude in Gjj minor from Book 1 of Das Wohltemperirte Klavier (BWV 863) in 
Ab minor. I do not believe this should be considered an error. What matters is whether 
or not the pitch interval names (e.g., 'rising major third', 'falling augmented fourth', etc.) 
between corresponding pairs of notes are the same in the computed spelling as they are in 
the published score. My implementation of Longuet-Higgins's algorithm therefore actually 
generates three alternative spellings for each piece: 1. the spelling S for the whole input 
passage computed by the algorithm; 2. the spelling that results when 3 is transposed by 
a rising diminished second; and 3. the spelling that results when S is transposed by a 
falling diminished second. Each of these three computed spellings is then compared with 
the pitch names in the original score and the number of errors made by the program is 
taken to be the number of errors in the best of the three alternative spellings. 

When my implementation of Longuet-Higgins's algorithm was run on the test "corpus 
described above, it made only 265 errors— that is, it predicted the correct pitch name for 
99.36% of the notes. 

Cambouropoulos's algorithm Unlike Longuet-Higgins, Cambouropoulos has not 
published an implementation of his pitch spelling algorithm, nor was he able to pro- 
vide me with his own implementation when I requested it. I therefore implemented my 
own version of his method, based on his published descriptions of it (Cambouropoulos, 
1996, 1998, 2000, 2001, 2002). 6 

Cambouropoulos's method involves first converting the input representation into a 
sequence of chromas or pitch classes in which the chromas are in the order in which 
they occur in the music (the chromas of notes that occur simultaneously being ordered 
arbitrarily). The chroma of a note can be computed from its MIDI note number using 

6 The LLsp code for my implementation of Cambouropp.ilos's pitch spelling algorithm Ls available online 
at nttp : //www . titanmusic . com/software . html. 
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the formula 

CHROMA = MIDI NOTE NUMBER mod 12. 

The term chroma has been used in this sense since around 1950 (Bachem, 1950; Burns 
and Ward, 1982, pp. 246, 262-264; Cross et al y 1991, pp. 212, 223-224; Deutsch, 1982, 
p. 272; Deutsch, 1999, p. 350; Shepard, 1964; Shepard, 1965; Shepard, 1982, p. 352; 
Dowling, 1991, p. 35; Ward and Burns, 1982, p. 432-433 ). The term chroma is essentially 
S3monymoiis with the term pitch class as used in a.tonal theory (Babbitt, 1960; Babbitt, 
1965; Forte, 1973; Morris, 1987; Rahn, 1980). 

Having derived an ordered set of chromas from the input, Cambouropoulos's algorithm 
then processes the music a window at a time, each window containing a fixed number of 
notes (set to a value between 9 and 15). Each window is positioned so that the first third 
of the window overlaps the last third of the previous window. Cambouropoulos allows 
'white note 5 chromas (i.e., 0, 2, 4, 5, 7, 9 and 11) to be spelt in three different ways 
(e.g., chroma 0 can be spelt as B}|, Ctj or Dbb) arid 'black note' chromas to be spelt in 
two different ways (e.g., chroma 6 can be spelt as F{j or Gb). Given these restricted sets 
of possible pitch names for each chroma, the algorithm computes all possible spellings 
for each window. A penalty score is then computed for each of these possible window 
spellings. The penalty score for a given window spelling is found by computing a penalty 
value for the interval between each pair of notes in the window and summing these penalty 
values. A given interval in a particular window spelling is penalised more heavily if it is 
an interval that occurs less frequently in the major and minor scales. An interval is also 
penalised if either of the pitch names forming the interval is . a double-sharp or a double- 
flat. For each window, the algorithm chooses the spelling that has the lowest penalty 
score. 

When I ran my implementation of Cambouropoulos's method on the test corpus, it 
made 2599 mistakes — that is, it predicted the correct pitch name for only 93.74% of the 
notes. When Cambouropoulos ran his own implementation of his method on the test 
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corpus, he found that it made even more errors than my implementation. This may be 
due to the fact that my implementation generates three alternative transpositions of the 
computed spelling and chooses the one that results in the least number of 



errors. 



Temperley's algorithm Temperley's (1997, 2001) pitch spelling algorithm is imple- 
mented in his harmony program which forms one component of his and Sleator's Melisma 
system. 7 The input to the harmony program must be in the form of a '^note-list" (Tem- 
perley, 2001, pp. 9-12) giving the MIDI note number of each note together with its onset 
time and duration in milliseconds. The harmony program also requires a specification 
of the metrical structure of the input passage which can be generated by first running 
the note-list through Temperley and Sleator's meter program (another component of the 
Melisma system). 

Temperley's (2001, pp. 115-136) pitch spelling algorithm searches for the spelling 
that best satisfies three "preference rules". The first of these rules stipulates that the 
algorithm should "prefer to label nearby events so that they are close together on the 
line of fifths" (Temperley, 2001, p. 125). This rule bears some resemblance to the basic 
principle underlying Longuet-Higgins's algorithm (see above). The second rule expresses 
the principle that if two tones are separated by a semitone and the first tone is distant from 
the key centre, then the interval between them should preferably be spelt as a diatonic 
semitone rather than a chromatic one (Temperley, 2001, p. 129). The third preference rule 
underlying Temperley's algorithm steers the algorithm towards spelling the notes so that 
a "good harmonic representation" results (Temperley, 2001, p. 131), a "good harmonic 
representation" being one that allows Temperley's harmony program to generate a correct 
harmonic analysis of the passage. 

Because the output of the harmony program depends on tempo, each OPND repre- 
sentation in the test corpus had to be supplemented with a file giving a tempo map for 
the piece represented . These two files were then automatically converted into a note-list 

7 The complete Melisma system together with documentation is available online at 
http : //www. link. cs . cmu.edu/music-analysis/. 
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in the form required by Temperley's programs. The output of the harmony program 
was then automatical^' converted into a format which would allow automatic comparison 
with the test corpus files. Again, for each piece in the test corpus, three transpositions of 
the spelling generated by Temperley's program were compared with the original spelling 
and the one with the least number of errors was selected. When Temperley's harmony 
program was run on the test corpus, it made only 122 mistakes — that is, it predicted the 
correct pitch name for 99.71% of the notes. 

Summary of the invention 

The invention described here consists of an algorithmic method called pslS that reliably 
computes the correct pitch names (e.g., Cjj4, Bb5 etc.) of the notes in a passage of tonal 
music, when given only the onset-time and MIDI note number of each note in the passage. 

The pslS algorithm has been shown to be more robust than previoxis algorithms, cor- 
rectly predicting the pitch names of 99.81% of the notes in a test corpus containing 41544 
notes and consisting of all the pieces in the first book of J. S. Bach's Das Wohltemperirte 
Klavier (i.e., psl3 incorrectly predicted the pitch names of only 81 notes in this test 
corpus) . 

Three previous algorithms (those of Cambouropoulos (1996, 1998, 2000, 2001, 2002), 
Longuet-Higgins (1976, 1987, 1993) and Temperley (1997, 2001)) were run on the same 
corpus of 41544 notes. On this corpus, Cambouropoulos's algorithm made 2599 mistakes, 
Longuet-Higgins's algorithm made 265 mistakes and Temperley's algorithm made 122 
mistakes. As psl3 made only 81 mistakes on the same corpus, this provides evidence in 
support of the claim that pslS is more robust than previous algorithms that attempt to 
perform the same task. 

The pslS algorithm is best understood to be in two parts, Part I and Part II. Part I 
consists of the following steps: 

1. computing for each pitch class 0 < c < 11 and each note n in the input, a pitch 
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letter name S(c,n) e {A, B, C, D, E, F,G}, calculated on the assumption that c is 
the tonic at the point in the piece where n occurs; 

2. computing for each note n in the input and each pitch class 0 < c < 11 a value 
CNT(c, n) giving the number of ti 

that includes n, some specified number /T pre of notes immediately preceding n and 
some specified number K post of notes immediately following n; 

3. computing for each note n and each letter name I, the set of chromas C(n,l) = 
{c.\S(r.,n)=l}- 

4. computing N(l,n) = ^ CNT(c, n) for each note n and each pitch letter name /■ 

ceC(n,l) 5 

5. computing for each note n, the letter name / for which N(l,n) is a maximum. 

Part II of the algorithm corrects those instances in the output of Part I where a 
neighbour note or passing note is erroneously predicted to have the same letter name as 
either the note preceding it or the note following it. Part II of pslS 

1. lowers the letter name of every Lower neighbour note for which the letter name 
predicted by Part I is the same as that of the preceding note; 

2. raises the letter name of every upper neighbour note for which the letter name 
predicted by Part I is the same as that of the preceding note; 

3. lowers the letter name of every descending passing note for which the letter name 
predicted by Part I is the same as that of the preceding note; 

4. raises the letter name of every descending passing note for which the letter name 
predicted by Part I is the same as that of the following note; 

5. lowers the letter name of every ascending passing note for which the letter name 
predicted by Part I is the same as that of the following note; 
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6. raises the letter name of every ascending passing note for which the letter name 
predicted by Part I is the same as that of the preceding note. 
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Detailed description of preferred implementations 



The invention consists of an algorithm, called pslS, that takes as input a representation of 
a musical passage (or work or set of works) in the form of a set / of triples (t,p c , d), each of 
"~ to " '" 6 ulK " UX1C,C0 L > Wlfc; cnrvmaitc pitch p c and the duration d of a single 
note or sequence of tied notes. The chromatic pitch of a note is an integer indicating the 
interval in semitones from AljO to the pitch of the note (i.e., Pc = MIDI NOTE NUMBER— 
21). The set / may be trivially derived from a MIDI file representation of the musical 
passage. 

If ei = (ti,p c ,i, d^ and e, = (tj,p Cij , dj) and e u e, € / then a is defined to be less than 
e,, denoted by e, < e,, if and only if U < t, or ft = tj A p Cji < PcJ ) or ft = tj A p c>i = 
Pcj A tfc < dj). Also, ei < ej - if and only if a < ej or a = ej . 

The first step in psl3 is to sort the set / to give an ordered set 



J~ ((tuPc,!^!) , (t 2 ,p Cr2 ,d 2 ) , . . . (t W , PcW ,d m )) (1) 

containing all and only the elements of /, sorted into increasing order so that j > i => 
(ti,Pa, di) < (tj,p CJ , dj) for all ft,p c>i , d,} , (t i>PBjl d,) e J. 
The second step in psi5 is to compute the ordered set 

C = <ci,c 2 ,...c fc ,...c|j|) ( 2 ) 

where c k = p cM mod 12, p Cjfc being the chromatic pitch of the fcth element of J. c k is the 
chroma of p ck . 

If A is an ordered set of elements, 

A — (a 1 ,a 2 ,... a k , . . .a\ A \) 
then let A[j] denote the (j + l)th element of A (e.g., A[0) = a u A[l] = a 2 ). Also, let 
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A\j,k] denote the ordered set that contains all the elements of A from A[j] to A[k — 1], 
inclusive (e.g., ^4[1.4] = {a2.a3.a4)). 

In addition to the set /, psl3 takes as input two mamerical parameters, which I call 
the precontext, denoted by i\ pre , and the postcontext, denoted by A post . The precontext 
must be an integer greater than or equal to 0 and the postcontext must be an integer 
greater than 0. The third step in pslS is to compute the ordered set 

$ = <*(C[0]), 0(C[1]), . . . 4>{C[k}), . . . 4>{C[\C\ - I))) (3) 

where 

4>{C[k]) = (CNT(0, k), CNT{\, fc), . . . CATT(11, fc)> (4) 

and CNT(c, k) returns the number of times that the chroma c occurs in the ordered set 
C[max({0, k — iT pre }),min({|C f |, k + impost})] (the functions max(B). and min(J3) return 
the greatest and least values, respectively, in the set B). 

Every pitch name has three parts: a letter name which must be a member of 
the set {A,B,C,D,E,F,G}; an inflection which must be a member of the infinite set 

f .... lim . U U L UU LLL LULL 1, J . _ , 1. ^ _ , %. • i, c i i. _ . • , . t-> 

. . pA, a, q, p 5 i/i/ , i/i/P, PPi/P, . . . j; cum an octave ituntuer wnicn must ue an integer, oy 
convention, if the inflection of a pitch name is equal to \\ it may be omitted, for example, 
Ct]4 may be written C4. The octave number of 'middle C 5 is 4 and the octave number of 
any other C is one greater than the next C below it on the staff and one less than the 
next C above it on the staff. The octave mamber of any pitch name whose letter name is 
not C is the same as that of the nearest C below it on the staff. Thus Bx3 sounds one 
semitone higher than C4 and Cb4 has the same sounding pitch as Bfcj3. 

If TV is a pitch name with letter name l(N) and octave number o(iV), then the morph 
of TV, denoted by m(iV), is given by the following table 



l(N) 


A B C D E F G 


m(N) 


0 1 2 3 4 5 6 
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The morphetic octave of N, denoted by o m (N), is given by 



o m (N) = 



o(N), if m = 0 orm = 1, 
o(N) — 1, otherwise. 



(5) 



The morphetic pitch of N, denoted by p m (N), is given by 

ftn(JV) = m(N) + 7o m (N). ( 6 ) 

The fourth step in pslS is to compute what I call the spelling table, S. This is done 
using the algorithm expressed in pseudo-code in Figure 4. In this pseudo-code, block 
structure is indicated by indentation and the symbol <W i s used to denote assignment 
of a value to a variable (i.e., the expression «x «- y» means "the value of variable x 
becomes equal to the value of y"). The symbol «©- is used to denote concatenation of_ 
ordered sets. Thus, if A and B are two ordered sets such that A = (a lt a 2 , ., . a ]A] ) and 
B = (61,62, • then 



/i ® if = 



(a x , 02, . . . am, £>!, 6 2 , . . . 6|b|) • 



Also, A®() = ()@A = A. 

Having computed the spelling table S, the algorithm goes on to compute what I call 
the relative morph list, R. This is computed using the algorithm in Figure 5. If A is an 
ordered set of ordered sets, then A\i]\j) denotes the (J + l)th element of the (i + l)th 
element of A. The expression S[k]\i] in line 6 of Figure 5 therefore denotes the (i + l)th 
element of the (k + l)th element of the spelling table S. If L is an ordered set and i is 
the value of at least one element of L, then the function POS(i,L) called in line 10 of 
Figure 5 returns the least value of k for which L[k] = i. 
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The next step in pslS is to compute the initial morph m mit which is given by 



rrhnit = <9init[C[0]] 



(7) 



where Q init = (0, 1, 1, 2, 2, 3, 4, 4, 5, 5, 6, 6). 

Next, psl3 computes the morph list, M, which satisfies the following condition 



(|M| = |C|) A (M[i] = (R[i] + m init ) mod 7 for all 0 < i < \C\). 



(8) 




(|0| = \C\) A (C>[i] - (J[i][0], C[t], M [i]> for all 0 < i < |C|). 



(9) 



Next, the ordered set ff , which I call the chord list is computed using the algorithm 
in Figure 6. 

In the next three steps, the algorithm ensures that neighbour note and passing note 
figures axe spelt correctly. Figure 7 illustrates the six types of error corrected by the 
algorithm. Figure 7a shows a lower neighbour note figure that is incorrectly spelt because 
the neighbour note (the middle note) has the same morphetic pitch as the two flanking 
notes. In this case, the morphetic pitch of the neighbour note is one greater than it should 
be. Figure 7b shows a similar error in the case of an upper neighbour note figure. In this 
case, the morphetic pitch of the neighbour note is one less than it should be. The passing 
note figures in Figure 7c and d are incorrect because the morphetic pitch of the passing 
note (the middle note) is one greater than it should be. Finally, the passing note figures 
in Figure 7e and f are incorrect because the morphetic pitch of the passing note is one 
less than it should be. 

psl3 first corrects errors like the ones in Figure 7a and b using the algorithm in 



Figure 8. If A is an ordered set of ordered sets then the expression ^4[z][7, k] denotes the 
ordered set (^[iJL?*], A[i]{j + 1], . . . A[i][k - 1]). Thus, the expression H[i + 2][fc][l, 3] in 
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line 4 of Figure 8 denotes the ordered set (H[i + 2][k][l], H[i + 2}[k}[2}). In Figure 7a and 
b, the neighbour note is one semitone away from the flanking notes. The algorithm in 
Figure 8 also corrects instances where the neighbour note is 2 semitones above or below 

the flanking notes but has the same morphetic pitch as the flanking notes. 

Next, pslS corrects errors like the ones in Figure 7c and e using the algorithm in 

Figure 9. Then it corrects errors like the ones in Figure 7d and f using the algorithm in 

Figure 10. In Figure 7c, d, e and f, the interval between the flanking notes is a minor third. 

The algorithms in Figures 9 and 10 also correct instances where the interval between the 



ilaiikmg notes is a major third. 



Having corrected neighbour note and passing note errors, pslS then computes a new 
morph list M' using the algorithm in Figure 11. 

Having computed M', it is now possible to compute a morphetic. pitch for each note. 
This can be done using the algorithm in Figure 12 which computes the ordered set of 
morphetic pitches, P. 

Finally, from J and P, pslS computes an OPND representation Z which satisfies the 
following condition 

(\Z\ = \J\) A (Z[i\ = (J\i][ 0 ], PPN({ J[i][l], />[*])), m2] ) for all 0 < i < \Z\). (10) 

The function PPN({ Pc , Pin )) returns the unique pitch name whose chromatic pitch is p c 
and whose morphetic pitch is p m . This function can be computed using the algorithm in 
Figure 13. In this algorithm, anything in double inverted commas (<<?") is a string— that 
is, an ordered set of characters. If A and B are strings such that A = «abcdef» and 
B = "ghijkr, then the concatenation of A onto B, denoted by A ®B, is "abcdefghijkF. 
In the pitch names generated by PPAT, the sharp signs are represented by V characters and 
the flat signs are represented by <f ' characters. A double sharp is represented by the string 
"ss". For example, the pitch name Cx4 is represented in PPN by the string "Css4". The 
empty string is denoted by "" and the function STR(n) called in line 23 returns a string 
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representation of the number n. For example, STR(-o) = " - 5", STR( 105. 6) = "105.6". 



Results of running pslS on the test corpus As explained above, pslS takes as input 
a set I of triples, (£,p c , d), each one giving the onset time, chromatic pitch and duration 
of each note. In addition, pslS requires two numerical parameters, the precontext, A pre , 
and the post context, A post . 

In order to determine the values of K pre and K post that give the best results, pslS 
was run on the test corpus 2500 times, each time using a different pair of values 
(ifp re , impost) chosen from the set i(K pre , K post ) j 1 < K pre , K poflt < 50}. For each pair of 
values (if pre, impost), the number of errors made by pslS on the test corpus was recorded. 

pslS made fewer than 122 mistakes (i.e., performed better than Temperley's algorithm) 
on the test corpus, for 2004 of the 2500 (K pre , K post ) pairs tested (i.e., 80.160% of the 
(K pre , K post ) pairs). 

pslS performed best on the test corpus when K pre was set to 33 and K post was set 
to either 23 or 25. With these parameter values, psl3 made only 81 errors on the test 
corpus — that is, it correctly predicted the pitch names of 99.805% of the notes in the test 
corpus. 

The mean number of errors made by pslS over all 2500 (K pre , K post ) pairs was 109.082 
(i.e., 99.737% of the notes were correctly spelt on average over all 2500 (K pTe , K post ) pairs). 
This average value is better than the result obtained by Temperley's algorithm for this 
test corpus. 

The worst result was obtained when both K pre and K post were set to 1. In this case, 
pslS made 1117 errors (97.311% correct). However, provided K pre is greater than about 
14 and K post is greater than about 21, pslS predicts the correct pitch name for over 
99.75% of the notes in the test corpus. 

A Lisp implementation of psl3 The Lisp implementation of psl3 given below as- 
sumes that the input representation / is represented as a list of sublists, each sublist 
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taking the form (t p c d) where t, p c and d are the onset time, chromatic pitch and du- 
ration, respectively, of a single note. For example, the passage in Figure 3a would be 
represented by the list 

((2 41 2) (4 43 2) (6 44 2) (8 46 2) (10 43 2) (12 44 1) (13 41 1) (14 40 1) 

^ v lu « <±, v2u <±o <±) ^24 48 5; (26 36 2) (28 38 2) (29 46 1) 

(30 39 2) (30 44 1) (31 43 1) (32 41 2) (32 46 1) (33 44 1) (34 38 2) 

(34 43 1) (35 41 1) (36 39 1) (36 43 2) (37 36 1) (38 35 1) (38 51 3) 

(39 36 1) (40 44 4) (41 50 1) (42 48 1) (43 50 1) (44 41 4) (44 50 1) 
(45 48 1) (46 47 1) (47 48 1)) 

Similarly, the output of the Lisp implementation of pslS given below is represented as 

a list of sublists. For example, the output of this imolerre-^— ~* 

7 u - laxo ^iwuwiuu,i.i V ii p^ io tor tne passage 

.given in Figure 3a is 

((2 »Dn4» 2) (4 »En4" 2) (6 »Fn4» 2) (8 "Gn4» 2) (10 »En4» 2) (12 »Fn4» 1) 

(13 »Dn4" 1) (14 »Cs4» 1) (15 »Dn4" 1) (16 "Bf4» 4) (20 »Gn4" 4) (24 »An4» 5) 

(26 "An3» 2) (28 »Bn3" 2) (29 »Gn4» 1) (30 »Cn4» 2) (30 »Fn4»' 1) (31 »En4" 1) 

(32 "Dn4» 2) (32 »Gn4» 1) (33 »Fn4» 1) (34 »Bn3» 2) (34 »En4" 1) (35 »Dn4" 1) 

(36 »Cn4» 1) (36 »En4» 2) (37 "An3» 1) (38 »Gs3» 1) (38 "CnS" 3) (39 »An3» 1) 

(40 «iFn4" 4) (41 "Bn4» 1) (42 »An4» 1) (43 »Bn4« 1) (44 »Dn4» 4) (44 »Bn4" 1) 

(45 "An4" 1) (46 "Gs4" 1) (47 "An4" 1)) 



Here, then, is the Lisp code for an implementation of pslS: 

(defun psl3 (^optional: (input -filename (choose-f ile-dialog) ) " 

(pre-context 33) 
(post-context 23)) 
(let* ( (sorted-input -representation 
(remove-duplicates 
(sort (with-open-file (input-file 

input-filename) 
(read input-file)) 
# ' vector-less-than) 
:test #>equalp)) 

(onset-list (mapcar #* first sorted-input-representation) ) 
. (chromatic -pitch-list 

(mapcar #' second sorted-input-representation)) 
(chroma-list 

(mapcar #' chromatic-pitch-chroma chromatic-pitch-list)) 
(n (list-length chroma-list)) 
( chroma- vector-list 

(do* ((cvl nil) 

(i 0 (1+ i))) 
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((= i n) 
cvl) 
(setf cvl 

(append cvl 

(list 
(do* ((context 

(subseq chroma-list 

(max 0 (- i pre-context) ) 
(min n (+ i post-context)))) 
(cv (list 00000000000 0)) 
(c 0 (+ 1 c))) 
((= c 12) 
cv) 

(setf (elt cv c) 

( count c cent ext )))))))) 
(chromamorph-table (list 0 1 122334556 6)) 
(spelling-table 
(do* ((f irst-morph nil nil) 
(spelling nil nil) 
(spelling2 nil nil) 
(st nil) 
(c 0 (1+ c))) 
((= c 12) 
st) 

(setf spelling 

(mapcar #' (lambda ( chroma- in-chroma-list) 
(elt chromamorph-table 

(mod (- chroma-in-chroma-list c) 12))) 
chroma-list) ) 
(setf f irst-morph (first spelling)) 
(setf spelling2 

(mapcar # ? (lambda (morph-in-spelling) 

(mod (- morph-in-spelling f irst-morph) 7)) 
spelling) ) 

(setf st (append st (list spelling2) ) ) ) ) 
(relative-morph-list 
(do ((morph- vector (list 0 0 0 0 0 0 0) 
(list 0 0 0 0 0 0 0)) 

(rml nil) 
(i 0 (1+ i)) 

(morphs-f or-this-chroma nil 

nil) 

) 

((= i n) 
rml) 

(setf morphs-f or-this-chroma 

(mapcar #' (lambda (spelling) 
(elt spelling i)) 
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spelling-table) ) 

(setf rml 

(do ((prev-score nil nil) 
(j 0 (1+ j))) 
((= j 12) 
; (pprint mo rph- vector) 
(append rml 

(list (position 

(apply #'max morph-vector) 
morph-vector) ) ) ) 

(setf prev-score 

(elt morph-vector 

(elt morphs-f or-this -chroma j))) 
(setf (elt morph-vector 



for this— clij. oma j)) 



(+ prev-score 

(elt (elt chroma- vector-list i) j))))))) 
(initial-morph (elt '(011223445566) 

(mod (first chromatic-pitch-list) 12))) 
(morph-list (mapcar #' (lambda (relative-morph) 

(mod (+ relative-morph initial-morph) 7)) 
relative-morph-list) ) 
(ocm (mapcar #'list onset-list chroma-list morph-list)) 
(ocm-chord-list (do* ((cl (list (list (first ocm)))) 

(i 1 (1+ i))) 
((= i n) 
cl) 

(if (= (first (elt ocm i)) 

(first (elt ocm (1- i)))) 
(setf (first (last cl)) 

(append (first (last cl)) 

(list (elt ocm i)))) 

(setf cl 

(append cl 

(list (list (elt ocm i)))))))) 
(number-of-chords (list-length ocm-chord-list)) 
; neighbour notes 
(ocm-chord-list 
(do* ((i 0 (1+ i))) 

((= i (- number-of-chords 2)) 
ocm-chord-list ) 
(dolist (notel (elt ocm-chord-list i)) 
(if (member (cdr notel) 

- (mapcar #'cdr (elt ocm-chord-list (+ i 2))) 
:test #'equalp) 
(dolist (note2 (elt ocm-chord-list (1+ i))) 
(if (= (third note2) 
(third notel)) 
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(progn 

(if (member (mod (- (second note2) (second notel)) 
12) 
J (l 2)) 
(setf (third note2) 

(mod (+ 1 (third note2)) 7))) 
(if (member (mod (- (second notel) (second note2)) 
12) 
'(1 2)) 
(setf (third note2) 

(mod (- (third note2) 1) 7)))))))))) 
; downward passing notes 
(ocm-chord-list 
(do* ((i 0 (1+ i))) 

( (= i (— number-cf— chords 2)) 
ocm-chord-list) 
(dolist (notel (elt ocm-chord-list i)) 

(dolist (note3 (elt ocm-chord-list (+ i 2))) 

(if (= (third note3) (mod (- (third notel) 2) 7)) 
(dolist (note2 (elt ocm-chord-list (1+ i))) 
(if (and (or (= (third note2) 
(third notel)) 
(= (third note2) 
(third note3))) 

« 0 

(mod (- (second notel) (second note2)) 
12) 

(mod (- (second notel) (second note3)) 
12))) 
(unless (remove-if 
#'null 
(map car 
#' (lambda (note) 

(/= (second note) 

(second note2))) 
(remove-if 
#'null 
(mapcar 
#> (lambda (note) 

(if (= (third note) 

(mod (- (third notel) 1) 
7)) 

note) ) 

(elt ocm-chord-list (1+ i)))))) 
(setf (third note2) 

(mod (- (third notel) 1) 7)))))))))) 

; upward passing notes 
(ocm-chord-list 
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(do* C(i o (1+ i))) 

((= i (- number-of -chords 2)) 
ocm-chord-list) 
(dolist (notel (elt ocm-chord-list i)) 

(dolist (note3 (elt ocm-chord-list (+ i 2))) 

(if (= (third note3) (mod (+ (third notel) 2) 7)) 
(dolist (note2 Celt 

ocm— chord— list (1+ i))) 
(if (and (or (= (third note2) 
(third notel)) 
(= (third note2) 
(third note3))) 

« 0 

(mod (- (second note3) (second note2)) 
12) 

(xod (- (second notc3) (second note!)) 
12))) 
(unless (remove-if 
#'null 
(mapcar 
#' (lambda (note) 

(/= (second note) 

(second note2))) 
(remove-if 
#'null 

(mapcar #> (lambda (note) 

(if (= (third note) 

(mod (+ (third notel) 
1) 

7)) 

note)) 

(elt ocm-chord-list (1+ i)))))) 
(setf (third note2) 

(mod (+ (third notel) 1) 7)))))))))) 
(morph-list (mapcar #>third (apply #>append ocm-chord-list))) 
(mor phe t i c -pi t ch- 1 is t 
(mapcar #' (lambda (chromatic-pitch morph) 

(let* ((morphetic-octavel (floor chromatic -pitch 12)) 
(morphetic-octave2 (+ 1 morphetic-octavel)) 
(morphetic-octave3 (- morphetic-octavel 1)) 
(mpl (+ morphetic-octavel (/ morph 7))) 
(mp2 (+ morphetic-octave2 (/ morph 7))) 
(mp3 (+ morphetic-octave3 (/ morph 7))) 
(chroma (mod chromatic -pitch 12)) 
(cp (+ morphetic-octavel (/ chroma 12))) 
(difference-list (list (abs (- cp mpl)) 

(abs (- cp mp2)) 
(abs (- cp mp3)))) 
(morphetic-octave-list (list morphetic-octavel 
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morphetic-octave2 
morphetic-octave3) ) 

(best-morphetic-octave 
(elt morphetic-octave-list 

(position (apply #'min difference-list) 
difference-list) ) ) ) 
(+ (* 7 be s t -mo rphetic -octave) morph))) 
chromatic-pit ch-list 
morph-list) ) 

(opd (mapcar # > (lambda (tpcd-triple morphetic-pitch) 
(list (first tpcd-triple) 

(list (second tpcd-triple) 

morphetic-pitch) 
(third tpcd-triple))) 



morphet i c-pi t ch-1 i st ) ) 
(opnd (mapcar #' (lambda (opd-datapoint) 

(list (first opd-datapoint) 

(p-pn (second opd-datapoint)) 
(third opd-datapoint))) 

opd))) 

opnd) ) 

(defun chromatic -pitch- chroma (chromatic-pitch) 
(mod chromatic-pitch 12)) 



(defun vector-less-than (vl v2) 
(cond ((null v2) nil) 
((null vl) t) 

((< (first vl) (first v2)) t) 

((> (first vl) (first v2)) nil) 

(t (vector-less-than (cdr vl) (cdr v2))))) 

(defun p-pn (p) 

(let* ((m (p-m p) ) 

(1 (elt '("A" M B" "C H "D" "E" "F" M G " ) m)) 
(gc (p-gc p)) 

(cdash (elt '(0235 78 10) m)) 

(e (- gc cdash)) 

(i 

(i (cond ((< e 0) 

(dotimes (j (- e) i) 

(setf i (concatenate 'string i "f")))) 
((> e 0) 
(dotimes (j e i) 

(setf i (concatenate 'string i "s' 1 )))) 
((= e 0) "n"))) 
(om (p-om p)) 



(oasa (if (or (= m 0) (= m 
om 

(+ 1 om))) 
(o (format nil "~D" oasa))) 
(concatenate 'string 1 i o))) 

(defun p-om (p) 
(div (p-pm p) 7)) 

(defun p-pm (p) 
(second p) ) 

(defun div (x y) 
(int (/ x y))) 

(defun int (x) 

(values (floor x))) 

(defun p-gc (p) 
(- (p-pc p) 

(* 12 (p-om p)))) 

(defun p-pc (p) 
(first p)) 

(defun p-m (p) 

(bmod (p-pm p) 7)) 

(defun bmod (x y) 
(- x 
(* y 

(int (/ x y))))) 
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CLAIMS 

1. A method for computing the pitch names (i.e., Cfj4, Bb3, etc.) of notes in a rep- 
resentation of music in which at least the onset time and MIDI note number (or 
chromatic pitch) of each note is given, comprising the steps of 

(a) computing for each pitch class 0 < c < 11 and each note n in the input, a pitch 
letter name S(c, n) G {A, B y C, D, E, F, <?}, calculated on the assumption that 
c is the tonic at the point in the piece where n occurs; 

(b) computing for each note n in the input and each pitch class 0 < c < 11 
a value CNT(c y n) giving the number of times that c occurs within a context 
surrounding n that includes n, some specified number K pTe of notes immediately 
preceding n and some specified number iv pos t of notes immediately following 
n- 

(c) computing for each note n and each letter name /, the set of chromas C(n, I) = 
{c\S(c,n)=l}; 

(d) computing N(l,n) = £ CNT(c,n) for each note n and each pitch letter 

cGC(n,l) 

name I , 

(e) computing for each note n, the letter name I for which N(l,n) is a maximum. 

2. The method of Claim 1, adapted to correct those errors in the output of the method 
of Claim 1 in which a neighbour note or passing note is erroneously predicted to 
have the same letter name as either the note preceding it or the note following it, 
comprising the further steps of 

(a) lowering the letter name of every lower neighbour note for which the letter 
name predicted b}' the method of Claim 1 is the same as that of the preceding 
note; 

(b) raising the letter name of every upper neighbour note for which the letter name 
predicted by the method of Claim 1 is the same as that of the preceding note; 
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(c) lowering the letter name of every descending passing note for which the letter 
name predicted by the method of Claim 1 is the same as that of the preceding 
note; 

(d) raising the letter name of every descending passing note for which the letter 
name predicted by the method of Claim 1 is the same as that of the following 
note; 

(e) lowering the letter name of every ascending passing note for which the letter 
name predicted by the method of Claim 1 is the same as that of the following 
note; 

(f) raising the letter name of every ascending passing note for which the letter 
name predicted by the method of Claim 1 is the same as that of the preceding 
note. 

Computer software or hardware adapted for performing the method of any preceding 
Claim 1-2. 
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ABSTRACT 

METHOD OF COMPUTING THE PITCH NAMES 
OF NOTES IN MIDI-LIKE MUSIC 
REPRESENTATIONS 

The invention described here consists of an algorithmic method called pslS that reli- 
ably computes the correct pitch names (e.g., CJj4, Bb5 etc.) of the notes in a passage of 

: ^ . — 1 . — i__ j.1 j. j_: ^ "* ,ttj~\ t j. _ i _ r . i_ y _ ;„ j.i. _ 

uwj.i<xi muoic, wiicii ^ivcii un^y unc oubct-biiijc tLiiu iviii^x liutt: u uniuti! ui efctcu no be ill Liie 

passage. 

The psi,? algorithm has been shown to be more robust than previous algorithms, cor- 
rectly predicting the pitch names of 99.81% of the notes in a test corpus containing 41544 
notes and consisting of all the pieces in the first book of J. S. Bach's Das Wohltemperirte 
Klavier (i.e., pslS incorrectly predicted the pitch names of only 81 notes in this test 
corpus). 

Three previous algorithms (those of Cambouropoulos (1996, 1998, 2000, 2001, 2002), 
Longuet-Kiggins (1976, 1987, 1993) and Temperley (1997, 2001)) were run on the same 
corpus of 41544 notes. On this corpus, Cambouropoulos's algorithm made 2599 mistakes, 
Longuet-Higgins's algorithm made 265 mistakes and Temperley's algorithm made 122 
mistakes. As pslS made only 81 mistakes on the same corpus, this provides evidence in 
support of the claim that psl3 is more robust than previous algorithms that attempt to 
perform the same task. 
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Figure 1: Examples of notes with identical MIDI note numbers being spelt differently in 
different tonal contexts (from Piston, 1978, p. 8). 
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Figure 2: Should the Ebs be spelt as Djis? (From Piston, 1978 ; p. 390.) 
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(a) 



(b) 



(2 D4 2), 
(13 D4 1), 
(29 G4 1), 
(35 D4 1>,, 
(44 B4 1), 
(30 C4 2), 
(39 A3 1), 



(4 E4 2), 
(14 CS4 1), 
(30 F4 1), 
(36 E4 2), 
(45 A4 1), 
(32 D4 2 ), 
(40 F4 4), 



(6 F4 2), 
(15 D4 1), 
(31- E4 l)v: , 
(38 CS3>, ' 
(46 GS4 1), 
(34 B3 2), 
(44 D4 4) 



■ y r r r 

(8 G4 2), 
(16 BF4 4), 
r ; - : ;(32 64 1)\ 
(41 B4 1), 
(47 A4 1>, 
(36 C4 1), 
} 



bet 

(10 E4 2>, 
(20 G4 4), 
:(33 F4 1), 
(42 A4 1), 
(26 A3 2), 
(37 A3 1), 



(12 F4 1), 
(24 A4 5), 
(34 E4 1), 
(43 B4 1), 
(28 B3 2>, 
(38 GS3 1), 



Figure 3: (a) Bars 1 to 4 of Bach's Fugue in D minor from Book 1 of Das Wohltempemrte 
Klavier (BWV 851). (b) The OPND representation of the score in (a). 
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1 si-o 

2 Q *- (0, 1, 1, 2, 2, 3, 3, 4, 5, 5, 6, 6) 

3 for c <— 0 to 11 

4 mo t- 0 

5 •>'>';" .;;*^>.-v.; : :;r r : vsl sac* 

6 *'<-<) 

7 for i<-0 to \C\ - 1 

8 5 <- s © (<2[(C[*"] - c) mod 12]) 

9 mo t— s[0] 

10 for z <- 0 to |C| - 1 

11 s' <- s' © ((s[i] - mo) mod 7) 

1A O O (±J (S ) 

13 return S 



Figure 4: Algorithm for computing the spelling table S. 
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1 *«-<) 

2 for % <-.0 .to |C| - 1 

3 V ^(6, 0,6, 0, 0, 0, 0) 
4 

5 for k <- 0 to 11 

7 for j <- 0 to 11 

8 5 prev <- V[/x[j}] 

10 R^-R@ (POS(max(K), K)) 

11 return it 



Figure 5: Algorithm for computing the relative morph list R. 
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1 H - «O[0]>> 

2 for i <37.1,.to: |C|-.l „' ; 

3 if b(i}[6}~= o\i - 1][0]' " 

4 tf[|#|-l]<-H[|tf|-l]e<0[*]> 

5 else 

6 If <- # e «0[i]» 

7 return H 



Figure 6: Algorithm for computing the chord list H. 
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$vf t r i i r i 'I 1 i r r I i !j ^ » fj i 

Figure 7: Neighbour note and passing note errors corrected by psi3. 
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1 for i^O to \H\ -3 



2 C-{) 

3 for k ^0 to \H[i + 2)\-l 

4 C-Ce;{i?^^][fc]fi;3]) - .- - 

5 for m ^'0 to -1 i 

6 if J7[tJ[ni][l,3] G C 

7 for n 2 <- 0 to + 1]| - 1 

8 if JT[i+l][n 2 ][2] = ff[i][n 1 ][2] 

9 if (£T[» + l][na][l] -iTMInJfl]) mod 12 € {1,2} 

10 H[i + l][n a ][2] 4- + l][n 2 ][2] + 1) mod 7 

11 if (J7[t][m][l] - H[i+ l][n 2 ][l]) mod 12 € {1,2} 

12 fl-[t + l][n 2 ][2] <- (fT[i + l][n 2 ][2] - 1) mod 7 



13 return H 



Figure 8: Algorithm for correcting neighbour note errors. 
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1 for i <- 0 to \H\ - 3 



2 for to \H[i]\ - 1 

3 for 7i 3 ^0 to \H[i + 2]| - 1 

4 if [i + 2][n 3 ][2] = (if \i][nx][2] - 2) mod 7 

5 for n 2 0 to + 1]| - 1 

6 if ir[t+l][fi2][2]=-J5r[ilfoiH?] or H[i -H l][n 2 ][2] = 2][n 3 ][2] 

7 if 0 < (H[i][m][l] ^ : Hli+ l}{n 2 ]{l]y mod 12 , < [n'i] [1] - H[i + 2][n 3 ][l]) mod 12 

9 for j <— 0 to - 1 

10 if + 1][?][2] = (iT[i][ni][2] - 1) mod 7 

11 C«-C©<J*[t + i]b"]> 

12 

13 for j ^ 0 to |C| - 1 

14 if ff[» + l][n 2 ][l]^Cb1[l] 

16 if 0 == () 

17 + i]N][2] <- (^WlmJM - 1) mod 7 



18 return i7 



Figure 9: Algorithm for correcting descending passing note errors. 
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1 for i <- 0 to \H\ - 3 



2 for ri! +- 0 to \H[i}\ - 1 

3 for n 3 «- 0 to |IT[t + 2] | - 1 

4 if #[» + 2][n 3 ][2] - (jyHImJp] + 2) mod 7 

5 for n 2 «- 0 to \H[i + 1]| - 1 

6 if H[i + 1] [n 2 ] [2] = [m] [2] or + [n 2 ] [2] = £T [t + 2] [n 3 ] [2] 

7 if 0 < (fT[< + 2][n 3 ][l] - H\i + l][n 2 ][l]) mod 12 <(lT[t + 2][n 3 ][l] - #MK][1]) mod 12 

8 C-0 • ■ ' ; ' ■ ' ' 

9 for j <-0 to |fT[t+'l]|-'i " 

10 if + [2] - (JIW[ni][2] + 1) mod 7 

11 <^Ce(H[i+l)[j]) 
12 

13 for j <- 0 to |C| - 1 

14 if ^[t + l][n 2 ][l]^C[j][l] 

is *«-*e<cb1> 

16 if 0= () 

17 H[i + l][n 2 ][2] <- (H[i][m]{2] + 1) mod 7 



18 return i7 



Figure 10: Algorithm for correcting ascending passing note errors. 
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1 o'-<> 

2 for i.<^ 0 -io \H\ -T ■ 

3 O' <- O'' 0 H[i] 

4 M' <— ( ) 

5 for i t- 0 to |0'| - 1 

6 M' <- M' ®{0'[i}[2]) 

7 return M' 



Figure 11: Algorithm for computing M' . 



THIS PAGE BLj^SK'Gjspto) 



12/13 



1 P«-<) 

2 for i <- 0 to | J\ - 1 

3 0l <- |J[i][l]/12j 

4 o 2 1 + Oi 

5 o 3 <— Ox — 1 

6 Pi <- Oj + M"'[i]/7 

7 p-2 <- o 2 + M'[i}/7 

8 p 3 «- o 3 + M'[i)/7 

9 c «- J[t][l] mod 12 

10 p'-e-o 1 +c/12 

11 £>-(|p'-Pii>'-p 2 |,|p'-p 3 |) 

12 a; «- (oi, 02, o 3 > 

14 P *- P © (M'[4 + 7o) 

15 return P 



Figure 12: Algorithm for computing P. 
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PPN(p) 

1 m <— p[l] mod 7 

2 L <— ( M" , "B" , "C" , "D" , "E" , "F" , "(?" ) 

3 I <- £[m] 

4 5c <-p[0] - 12(Lp[i]/7j) 

5 ,4 <- (0,2,3,5,7,8,10) 

6 c' -e- A[m] 

7 e g c - d 

8 i <- "" 

9 if e < 0 

10 for j f— 0 to -e - 1 

11 i ^i® "/" 

12 else 

13 if e > 0 ' = " ■ ' ■ 

14 for j <— 0 to e - 1 

15 i ^- i ® V 

16 else 

17 „' , «„" 

a. » t> * /t 

18 o m <- Lp[l]/7J 

19 if ra = 0 or ra — 1 

20 o<- o m 

21 else 

22 o ^ 1 + o m 

23 o str <- STH(o) 

24 return I 0 z 0 o 



Figure 13: PPiV algorithm. 
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